Classical Loschmidt echo in chaotic many-body systems 



in 
o 
o 

(N 



(N 

Q 

U 
d 



> 
m 

o 
m 
o 
in 
o 



> 

x 



Gregor Veble 1 ' 2 and Tomaz Prosen 1 
1 Physics Department, FMF, University of Ljubljana, Ljubljana, Slovenia 
2 Center for Applied Mathematics and Theoretical Physics, University of Maribor, Maribor, Slovenia 

(Dated: February 8, 2008) 

General theoretic approach to classical Loschmidt echoes in chaotic systems with many degrees 
of freedom is developed. For perturbations which affect essentially all degrees of freedom we find a 
doubly exponential decay with the rate determined by the largest Lyapunov exponent. The scaling 
of the decay rate on the perturbation strength depends on whether the initial phase-space density 
is continuous or not. 
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Understanding macroscopic irreversibility from re- 
versible miscroscopic evolution laws is one of the fun- 
dametal controversies of statistical mechanics, put for- 
ward by Loschmidt Q. The problem has been reformu- 
lated in terms of the so called Loschmidt echo or fidelity 
QilHOLIIlEli mainly in the context of quantum chaos and 
quantum information, and only more recently 0, H @ m 
the context of classical dynamics. In @ it was shown 
that dynamics in the echo picture, which is a compo- 
sition of perturbed forward and unperturbed backward 
Hamiltonian evolution in time, can again be written as 
forward Hamiltonian evolution. The Hamiltonian in this 
case becomes time dependent and is merely the perturba- 
tion part of the original Hamiltonian evaluated at a point 
obtained by the forward unperturbed evolution. Analyz- 
ing the motion in the echo picture for the case of Anosov 
systems showed that, for sufficiently long times, it can 
be separated into a composition of independent expo- 
nentially accelerated one-dimensional motions, each such 
motion locally corresponding to a Lyapunov vector direc- 
tion of the original, unperturbed system. Motion along 
each such (unstable) direction could be well described 
statistically, using a kernel whose shape is time indepen- 
dent but whose width increases as ~ e Ai * where A^ is the 
corresponding Lyapunov exponent. The theory shows 
that such an evolution leads to the Lyapunov decay of 
fidelity for the case of a system with one degree of free- 
dom, confirming results of Refs. 0, Q, @ • For systems 
with more than a single degree of freedom a cascade of 
exponential decays with increasing decay rates was pre- 
dicted and observed. 

Here we extend the results to systems with many de- 
grees of freedom. While a similar theroy may still be ap- 
plied, the results are significantly different as the regime 
which is considered perturbative in systems with few de- 
grees of freedom is actually the dominant one for systems 
with a large number of freedoms. What we find is that 
the decay of fidelity for systems with a large number of 
degrees of freedom decays as a doubly exponential func- 
tion of time for the case of a global perturbation (afftect- 
ing all degrees of freedom), and the scaling of decay rate 
on the strength of perturbation depends on whether the 
initial phase space density is continuous or not. In an- 
other case of a local perturbation we find and justify a 



system dependent form of decay that falls between the 
exponential and doubly exponential behaviors. 

In a system with N phase space variables fidelity is 
defined as the overlap 
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d x po(x) p t {x), 



where pg is the initial phase space density, a; is a phase 
space point and pt is the echo density as obtained by a 
composition of a perturbed forward evolution for time t 
with Hamiltonian H — Ho + SV, 5 is the strength of per- 
turbation, and an unperturbed backward evolution, again 
for time t, with Hamiltonian Hq. As shown in [ij, un- 
der the assumption that the time is long enough in order 
for the echo dynamics to start to behave statistically, the 
echo density can be written as 
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where new variables j/j are introduced (see 0) that give 
the coordinates along the stable and unstable local Lya- 
punov vectors. The kernels Ki are time independent 
with their typical widths of the order of 1 , the time de- 
pendence being present in the widths o~i = 5^% exp(A;t) 
for the unstable and cr, = 5"fi for the stable directions, 
where S is the strength of the perturbation and Aj the 
corresponding Lyapunov exponent. The factor 7, in the 
expression of Gi is determined by the typical values of 
the phase space variable Wi(x) = ei(x) ■ W(x), where 
e_j is the corresponding Lyapunov vector defined as in 
(3. Fidelity can be rewritten in the new coordinates as 
f(t) = Jd N y Po (y) Pt (y). 

Let us assume now that the initial density is a prod- 
uct po(y) — YiiLiPu 1 f° r example being a char- 
acteristic funciton on a hypercube, or a multidimen- 
sional Gaussian, in coordinates yi. In this case the 
subsequent analysis is very much simplified and we see 
no reason for the results to be different in the general 
cases, as indeed the numerics, which does not corre- 
spond to initial product denisty, shows. Under this as- 
sumption fidelity itself can be written as a product of 
contributions from different individual Lyapunov direc- 
tions, f(t) — Yli=i /»'(*)■ Each such contribution is of 
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the form = J dy p^ ] {y) p[ l) (y), where p t W (y) = 

j: J dy'p^ {y + y')K l (^j . The initial density contribu- 
tions po can be considered to be square normalized, and 
we may write 



Mt) = l~ I dydy' pf{y) 
La, 
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we 



Introducing Di(y') = J dy 

may see that the behavior of f(t) will be markedly differ- 
ent depending on whether the function po is continuous 
or not. In the continuous case the function Di(y') oc y' 2 
for small y', whereas in the discontinuous but piecewise 
continuous case the behavior is Di(y') — G, \y'\, where 
G represents the sum of all squares of the discontinuity 
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From this it follows that to the lowest order 
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(1) 



where [3 — 1 for the discontinuous case, and [3 = 2 for 
the continuous case, and Ci is a constant. 

In the case of many degrees of freedom fidelity is a 
product of contributions of many separate unstable di- 
rections. Even if individual contributions fi(t) do not 
deviate significantly from 1 that may not be the case for 
the product of all /,. When N becomes very large we may 
therefore expect the regime where equation is valid 
to describe the behavior of fidelity for all times until it 
becomes insignificantly small, at least in the cases where 
a significant fraction of unstable directions contribute to 
the decay of fidelity. The resulting expression is thus 
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(Sjif exp(J3Xit) 
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where the product is taken along all the unstable direc- 
tions i, Ci are positive coefficients depending on the shape 
of the initial density pa and the details of the perturba- 
tion V and (3 — 1,2, depending whether the initial den- 
sity is discontinuous, or continuous, respectively. 

The behavior of fidelity depends on a few consider- 
ations. Here we limit ourselves to the cases of initial 
phase space densities that are not constant (inpedepen- 
dent) with respect to any of the phase space variables. 
Another important distinction is the nature of the per- 
turbation. The perturbation can be either local in na- 
ture, meaning that it significantly affects only a finite 
number of degrees of freedom of the system, or they can 
be extended such that each degree of freedom is affected 
equally on average. 

In the case of an extended perturbation, each compo- 
nent of the gradient of the perturbation potential VV 
is of the order of 1, i.e. independent of N. As it may 
be argued that the norm of the Lyapunov vector is 
independent of N, the magnitude of the variable Wj 

can be estimated as (e$ • VV) 2 = fX^jLi [e»]^ djV 




0.4 0.6 

(i-1/2)/L 



FIG. 1: The Lyapunov spectrum of the kicked Ising chain as 
calculated after 10 4 kicks for L — 10 (full circles), 100 (empty 
circles) and 1000 (line) spins sorted from the largest to lowest 
positive value (the negative and values of the spectrum are 
omitted). On the horizontal axis we give (i — 1/2) /L, where 
i is the consecutive index of the Lyapunov exponent. 



Ej=i [ e *lj (®jV) 2 — ^(1)- The above estimate is done 
by assuming that the components of the Lyapunov field 
vector are not correlated. This estimate gives the mag- 
nitude of the factor 73 in equation (J2J. In that equation, 
the most relevant contributions to the decay are those 
with the highest Lyapunov exponents. As numerical evi- 
dence below suggests, for uniform systems the Lyapunov 
spectrum in the thermodynamics limit N — > 00 has a 
limiting distribution of Lyapunov exponents, with an ex- 
isting maximum exponent A max . From this it follows that 
in an interval [A max — AA, A max ] there will exist oc AXN 
Lyapunov exponents in this interval. The number n of 
relevant factors that contribute to the decay of fidelity as 
given in J2J therefore scales linearly with TV. Assuming 
that the factors c, (#7i)^ in equation have a well de- 
fined average and are statistically independent, when N 
becomes very large equation @ can be written as 



f(t) = exp (-aNS f3 exp(/?A max i)) 



(3) 



where j3 = 1 or 2. The decay of fidelity in many body 
systems with global perturbations is therefore expected 
to decay as a doubly exponential function, with the decay 
being stronger with increasing N. This is the main result 
of this Letter. 

When considering the case of a local perturbation more 
care needs to be taken. Because the perturbation is local, 
the variables Wj that determine the spreading of densi- 
ties along individual Lyapunov directions are highly cor- 
related amongst themselves. Namely, the Lyapunov vec- 
tors corresponding to the exponents close to the maxi- 
mum one have a very similar structure on the local scale 
(in the sense of neighbouring degrees of freedom) and 
therefore the corresponding different phase space vari- 
ables Wi are, according to their definition, all very sim- 
ilar. This means that, for very short times, the contri- 
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FIG. 2: Decay of fidelity as a function of time for the case 
of the discontinuous distribution. We divide the logarithm 
of fidelity by L, show it in log-scale, and shift the time 
by log(5)/A max . We show the data for L — 10 (triangles), 
100 (circles) and 1000 (squares) with perturbation strengths 
8 = 10~ 3 (empty symbols) and 10~ 5 (full symbols). The num- 
ber of trajectories was nt = 10 7 for the cases of L — 10, 100 
and nt = 10 6 for L = 1000. We may see that in the regime 
of validity of equation as shown by the dotted line for 
j3 = 1 with an arbitrarily chosen constant a, all data scale 
to the same curve in the diagram. The deviations near the 
saturation values are due to different times of dynamics of 
the system and therefore a different decay of correlations (see 
§), this regime however moves to the values of fidelity ap- 
proaching with L — ► oo. 



butions to the spreading of density along different Lya- 
punov directions cannot be considered to be independent 
of each other, but rather the spreading is effectively one- 
dimensional along the average of the Lyapunov directions 
in phase space. This average direction corresponds to the 
(local) degree of freedom in which the perturbation is be- 
ing applied. At longer times, due to the coupling between 
the neighbouring degrees of freedom, the differences be- 
tween different Lyapunov vectors and exponents become 
significant and the spreading of densities takes place in 
an increasing number of dimensions. This can be under- 
stood as the diffusion of the effect of the perturbation 
to the neighbouring degrees of freedom of the system. 
Thus we may conclude that: (i) the decay of fidelity for 
a local perturbation is independent of the size N of the 
system for sufficiently large systems, but (ii) that it is 
faster than exponential since more and more degrees of 
freedom are beginning to contribute to the decay due to 
the diffusion of the perturbation effects (exponential de- 
cay is expected if only a single degree of freedom would 
be relevant). However, since the nature of the diffusion 
and the related properties of the Lyapunov vectors and 
the spectrum are very much system dependent, there is 
no universal decay of fidelity for the case of the local 
perturbation. 

Our theoretical predictions will be demonstrated nu- 
merically in a generic chaotic classical many-body sys- 
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FIG. 3: Decay of fidelity as a function of time for the case 
of the continuous distribution, using the same represantation 
of scaled data as in Fig. H We show the data for L = 10 
(triangles) and 30 (circles) with perturbation strengths 8 = 
10~~ 3 (empty symbols) and 10 -4 (full symbols). The number 
of trajectories was n t = 10 7 . We show the prediction with 
arbitrarily chosen a and /3 = 2 using the dotted line. 



tern, namely a classical kicked Ising spin chain with the 
Hamiltonian (studied in the quantum version in p|) 



T q Z <J Z 



L classical spin variables s*j take values on unit spheres, 
sf = 1, and periodic boundary conditions are assumed 
L + 1 = 1. The coupling and kick parameters were nor- 
mally set to constant values Ji — 2 and K t — 2, respec- 
tively. The dimension of the phase space of such a system 
is N = 2L. The perturbation of the system was done by 
modifying the kicking strength K[ = Ki + SVi. In the 
case of the global perturbation this was done for all i 
uniformly, normally we set Vj = 1, while in the case of 
the local perturbation, Vi = 1, = 0. 

In order to demonstrate the differences between con- 
tinuous and discontinuous distributions, we used two sets 
of initial densities. For the discontinuous case /3 = 1, we 
used an initial uniform distribution of all spins on the 
upper hemispheres s z > 0. Even though the support for 
these distributions is not small in length as compared to 
a diameter of phase space, the fraction of the phase space 
volume that they cover, which gives the fidelity plateau 
/(oo), drops to as 2~ L in the thermodynamic limit. 

For the continuous density our choice of the distri- 
bution was a product of square normalized distribu- 
tions w(,s z ) = ^/3/8 (s z + 1) for all spins. This, how- 
ever, does not imply that the distribution is a prod- 
uct distribution in the Lyapunov directions. Similarly 
to the case of discontinuous densities, fidelity plateau is 
again expected to drop to in the large L limit, namely 
f(t -> oo) = [2- 1 j\ds z w{s*)} L = (3/8) L / 2 . 
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FIG. 4: Decay of fidelity as a function of time when a single 
spin is perturbed, and discontinuous initial density is cho- 
sen. We show the decay in the case of L = 100 (empty, 
nt ~ 5.8 x 10 6 trajectories) and L — 1000 (full triangles, 



n t = 10 b ) with S 
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The prediction for extended 



perturbations with arbitrary a is shown as a dotted line and 
the (arbitrarily shifted) exponential decay with the maximum 
Lyapunov exponent of the system as a dashed line. As ex- 
plained in the text, the actual decay falls in between the two 
types of behavior. 



To be able to test the predictions we need to deter- 
mine the Lyapunov spectrum of such a system. We show 
the positive part of the Lyapunov spectrum in Fig. ^ 
We checked that the negative part of the spectrum cor- 
responds to the positive part and that the spectrum con- 
verges with time. We may see that the actual density of 
Lyapunov exponents divided by L beautifully scales to 
the limiting density when L is large enough. Such a be- 
haviour is typically found for interacting many-particle 
systems, such as e.g. Fermi- Past a- Ulam chain, Lorentz 
gas, etc. This justifies our assumption that the number of 
Lyapunov exponents in the neighbourhood of the largest 
Lyapunov exponent increases linearly with N = 2L, yet 
the maximum Lyapunov exponent converges with N. 

We show the decay of fidelity for the discontinuous 
case and a global perturbation in Fig. [3 Using the 
appropriate scaling with respect to the number of spins 
L = N/2 and perturbation strength S as predicted in 



equation ©, where (3 — 1, we may see that the doubly 
exponential behavior with the predicted dependence on 5 
and N holds true. It should be noted that, while there is 
an observed plateau of fidelity when N is still small, this 
plateau exponentially tends to with N and indeed the 
doubly exponential behavior of fidelity can be observed 
for all values of fidelity for large N. 

Similar observations hold true for the continuous case 
where (3 = 2 as shown in Fig. |3| For the same system, us- 
ing continuous densities results in a faster form of decay 
with time and with quadratic dependence on the pertur- 
bation. In the continuous case we also observe expected 
numerical statistical fluctuations of fidelity at values close 
to 1, especially for large N. This is because the initial 
density is represented by a finite number of trajectories 
and even when t = the estimation of fidelity will differ 
from 1. Nevertheless, the predicted scaling with N and 
S is observed along with the faster decay of fidelity with 
time as compared to discontinuous case. 

For the case of a local perturbation with discontinuous 
densities as presented in Fig. 0] we observe a decay that 
lies in between the exponential behavior, as expected in 
systems with a single degree of freedom, and the doubly 
exponential behavior as is found in the global perturba- 
tion case. Details of such behavior are system dependent 
as they are determined by the diffusion of the effects of 
the perturbation to neighbouring degrees of freedom. 

In addition we have tested our predictions for the dis- 
ordered cases where the kicking strengths Ki randomly 
varied from spin to spin, as well as the perturbations of 
the kicking strength V%, however the results were again 
confirming our theory. 

In conclusion, we have shown and explained that the 
decay of the classical Loschmidt echo in systems with a 
large number of degrees of freedom decays as a doubly 
exponential function for a global perturbation, and the 
decay rate depends on whether the phase space densities 
chosen are continuous or not. For the case of a local per- 
turbation a system dependent form of decay is observed 
and justified that falls between exponential and doubly 
exponential behaviors. We acknowledge support by the 
grant P 1-044 of Ministry of Science and Technology of 
Slovenia, and in part by the U.S. ARO grant DAAD19- 
02-1-0086. 
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